ac Stark shift and multiphoton-like resonances in low-frequency driven optical lattices 
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We suggest that Bose-Einstein condensates in optical lattices subjected to ac forcing with a 
smooth envelope may provide detailed experimental access to multiphoton-like transitions between 
ac-Stark-shifted Bloch bands. Such transitions correspond to resonances described theoretically by 
avoided quasienergy crossings. We show that the width of such anticrossings can be inferred from 
measurements involving asymmetric pulses. We also introduce a pulse tracking strategy for locating 
the particular driving amplitudes for which resonances occur. Our numerical calculations refer to a 
currently existing experimental set-up [Haller et al., PRL 104, 200403 (2010)]. 

PACS numbers: 03.75.Kk, 42.50.Hz, 67.85.Hj 
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I. INTRODUCTION 

The study of multiphoton excitation and ionization 
processes so far has concerned the interaction of mat- 
ter with strong electromagnetic fields |l|-|4j. Extending 
our previous works 0, Q , we demonstrate in the present 
proposal that ideas and concepts developed in this field 
of research can also be applied for understanding the re- 
sponse of ultracold macroscopic matter waves in opti- 
cal lattices to ac forcing with a slowly varying envelope. 
Such systems offer additional handles of control which 
are not available in experiments with atoms or molecules 
exposed to pulses of laser radiation. Therefore, they can 
give novel insights into multiphoton dynamics in general, 
and in particular may also allow one to systematically in- 
vestigate the effects of interparticle interactions on such 
transitions. 

The experimental investigation of Bose-Einstein con- 
densates (BECs) in time-periodically forced optical lat- 
tices has been pushed forward with remarkable vigor 
within the last years, addressing quite diverse topics 
such as parametric amplification of matter waves Q ) dy- 
namic localization [8, 9], photon-assisted tunneling [Toj . 
coherent control of the superfluid-to-Mott insulator tran- 
sition fll |. quantum ratchets [12j, super Bloch oscil- 
lations I13I1 . quantum simulation of frustrated mag- 
netism [l4j , controlled correlated tunneling [la . 1 161 ] , and 
even the realization of tunable gauge potentials [l7j . 
Thus, BECs subjected to time-periodic forcing constitute 
an emerging branch of research [l8j . Here we outline how 
systematic use of the forcing's envelope will allow one to 
exctract crucial information about the systems' dynam- 



II. DETECTION OF MULTIPHOTON-LIKE 
RESONANCES 

The present theoretical study of multiphoton-like con- 
densate dynamics refers to conditions recently realized 
by Haller et al. [l3[: These authors have loaded BECs 
of Cs atoms into a vertically oriented ID optical lattice 
V{x) = (Vq/2) cos{2ki^x). Here k\, = 2tt/\l, where Al = 



1064.49 nm is the wavelength of the lattice-generating 
laser light, so that the lattice period is d = ir/k^. 
We choose a comparatively shallow lattice with depth 
V = 2.3 E T , where E r = (hk^) 2 / '(2m) is the single- 
photon recoil energy [l9l ]. with m denoting the mass of 
the atoms. The first excited band then is of the "above- 
barrier" type, so that atoms in this and all higher bands 
do not need to tunnel through the barriers when moving 
over the lattice. Using a magnetically induced Feshbach 
resonance, the s-wave scattering length of the Cs atoms 
is tuned to zero [2(| , so that one is dealing with a conden- 
sate of effectively noninteracting particles. By means of 
a time-periodic modulation of the levitation gradient em- 
ployed for compensating gravity, an external oscillating 
force is introduced which acts on the atoms with max- 
imum amplitude and frequency v = lj/(2tt), such 
that their dynamics are governed by the Hamiltonian 
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V(x) - s(t)F max xcos(ojt) 



(1) 



with a dimensionless shape function s(t) determining 
the envelope of the pulses applied. All calculations re- 
ported here are performed for the relatively low frequency 
w/(2tt) = 300 Hz, so that hw = 0.23 E T . As indicated in 
Fig.[TJ this means that the gap Ai = 1.14 E r between the 
lowest two Bloch bands Ei{k) and E^{k) of the unper- 
turbed optical lattice, which occurs at the Brillouin zone 
edge k/kh = ±1, amounts to 5.05 hoj: Exciting particles 
from the initially occupied lowest band to the first excited 
one requires the absorption of more than five "photons" . 

A first glimpse at the condensate dynamics under the 
action of a force F(t) is provided by Bloch's acceleration 
theorem: The expectation value (k)(t) of an atomic wave 
packet in fc-space evolves in time according to 



h-(k)(t) = F(t) 



(2) 



provided interband transitions remain negligible [211 ] . As- 
suming a sinusoidal force F(t) — F cos(wf)0(t) instanta- 
neously switched on at time t = 0, this gives (k)(t)/k-^ = 
(k)(0)/k L + (K/n)sm(u>t), where K = Fd/(Hu>) is a di- 
mensionless measure of the driving amplitude. Thus, 
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FIG. 1: Energy band structure of an unperturbed optical 
cosine lattice with depth Vo = 2.3 E Y . Measured in terms 
of the "photon" energy fiu = 0.23 E r , as corresponding to 
the driving frequency uj/(2n) — 300 Hz, the gap between the 
lowest two bands is Ai = 5.05 hw at the Brillouin zone edge, 
and Ao = 18.24 hw at its center. 
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FIG. 2: Interband transition probabilities after pulses with 
the squared-sine envelope ((3J and length T p = 60 x 2-k/oj, ob- 
tained from numerical solutions of the Schrodinger equation. 
All atoms which are excited to bands n > 1 after a pulse are 
assumed to escape from the lattice. 



Hamiltonian ([TJ, with a smooth envelope 

s(t) = sin 2 (7rf/T p ) ; < t < T p . (3) 

The pulse length is fixed at 60 driving cycles, T p = 60 T 
with T = 2-k/uj. In order to model the dynamics of a 
noninteracting BEC responding to such pulses, we start 
with an initial state 

1>(x, 0) = ^ J dk 9l (k, 0)xi,kfr) ( 4 ) 

made up from Bloch waves Xi,k(%) of the lowest band, 
employing a Gaussian momentum distribution 

9 ^,0)^(^Aky^e X ^JJ^^"j (5) 

centered around (fc)(0)/fcL = with width Ak/k^ = 0.1, 
and then solve the single-particle Schrodinger equation 
by means of a Crank-Nicolson algorithm [23|]. Vary- 
ing the maximum scaled amplitude K ma _ x = F max d /(Hoj) 
from pulse to pulse, we plot the escape probability from 
the lattice at the end of each pulse, at t = T p . Here 
we assume that only atoms which finally still populate 
the lowest band remain in the lattice, since atoms which 
have been excited to the higher "above-barrier" bands 
tend to escape from the shallow lattice quite fast. Ev- 
idently, interband transitions start to make themselves 
felt at K max m 2.5, and reach substantial strength 
when K max w 3, confirming the above rough estimate. 
Thus, proof-of-principle experiments performed along 
these lines should establish the feasibility of using BECs 
in driven optical lattices as novel probes for multiphoton- 
like transitions: Take driving pulses with smooth en- 
velopes, and measure the interband transition probability 
at the end of each pulse. In a series of such measurements 
with a constant pulse shape one then should observe a 
pronounced onset of interband transitions when the max- 
imum amplitude is successively increased. 



III. AVOIDED-QUASIENERGY-CROSSING 
SPECTROSCOPY WITH ASYMMETRIC PULSES 



when considering a wave packet initially at rest in the 
lowest band, so that (fc)(0)/fcL = 0, the packet's center 
(k)(t) just reaches the Brillouin zone edge when K = n. 
Because the zone edge gives rise to Zener-type transi- 
tions ■22j, one then expects pronounced excitation of 
higher bands. Therefore, one has K 7r as a rough 
order-of-magnitude estimate of the amplitude required 
for inducing multiphoton-like transitions. 

This expectation is confirmed by the numerical calcu- 
lations summarized in Fig. [5J Here we consider driving 
pulses s(t)F max cos(wi) as already incorporated into the 



In a second step, this approach can be employed for- 
getting more detailed insight into multiphoton dynam- 
ics. Namely, the single-band acceleration theorem @ 
ignores an essential element: Not only does the wave 
packet's center (k)(t) move within its band in response 
to the external forcing, but also the bands themselves 
are "dressed" by the drive, and therefore experience an 
ac Stark shift [2J]. Hence, initially nonresonant bands 
may be shifted such that their separation in energy ap- 
proaches an integer multiple of Huj for certain Bloch 
wavenumbers k, possibly leading to strong resonant in- 
terband coupling, that is, to a multiphoton resonance. 
Theoretically, such resonances are found by computing 
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FIG. 3: Quasienergies e„(0) at the center of the Brillouin zone 
for an optical lattice with depth Vo/E r = 2.3, driven with 
scaled frequency fku/E r = 0.23. The quasienergy originating 
from the lowest band n = 1 is marked by arrows. It exhibits a 
substantial ac Stark shift, and undergoes pronounced avoided 
crossings for K > 3. 




FIG. 4: Quasienergies e„(k) for k/ki, — 0.8 and n = 1,2,3. 
The arrow indicates the quasienergy originating from the low- 
est energy band n = 1. 

the quasienergy bands e n (k) which emerge from the en- 
ergy bands E n (k) in the presence of a drive with constant 
amplitude. These quasienergy bands reflect the ac-Stark- 
shifted energy bands, projected into an interval of width 
Ae = hu, so that a multiphoton resonance translates 
into an avoided quasienergy crossing [D-[6[ . For exam- 
ple, Fig.[3] shows the quasienergies e n (0) with n = 1,2, 3, 
which pertain to the pulses considered in Fig. [2] plotted 
versus the scaled amplitude K. The quasienergy originat- 
ing from the ground-state energy E\ (0) of the optical lat- 
tice undergoes two well-resolved avoided crossings when 
K > 3, signaling the presence of two individual multi- 
photon resonances. The observation that these "large" 
resonances begin to show up only for K « it nicely re- 
lates the elaborate quasienergy approach to the previous 
elementary reasoning based on Eq. @. 

Inspecting that elementary reasoning once again, one 



expects multiphoton resonances to occur for smaller driv- 
ing amplitudes when the initial packet is centered around 
a nonzero wavenumber, (fc)(0)/fcL 7^ 0, since then smaller 
amplitudes are required for reaching the Brillouin zone 
edge. This expectation is confirmed by Fig. [4j which 
depicts quasienergies e„(fc) for kjk\, = 0.8, with the 
first resonance showing up already at K w 0.9. Ex- 
perimentally, one can prepare an initial state with ar- 
bitrary (fc) (0) by subjecting the condensate to a suitable 
"kick" |25|, |2fj. Thus, one should also be able to detect 
the resonances predicted by Fig. 2J To this end, we pro- 
pose a particular kind of "avoided-quasienergy-crossing 
spectroscopy" based on the use of asymmetric pulses sit). 
For illustration, we assume that the rising part of such 
pulses be given by the first half of the envelope ((3|) , with 
fixed switch-on time Tp 1 '/2 = 5T, while their decreasing 

part is described by the second half of a squared-sine en- 

(2) 

velope, but with a different switch-off time Tp /2. Let 
the maximum scaled amplitude be K mayi — 1.2. During 
the rising part of such a pulse, a wave packet initially cen- 
tered around (h)(0) /ki, = 0.8 then follows its quasienergy 
states adiabatically, until the instantaneous amplitude 
reaches the multiphoton resonance at K ss 0.9 visible 
in Fig. |4j Then the packet undergoes a Landau-Zener 
transition to the anticrossing quasienergy state [27]. Due 
to the rapid switch-on of the pulse, and to the narrow 
quasienergy separation Se at the avoided crossing, that 
transition is almost complete. Thereafter, the packet 
again adiabatically follows the pulse envelope, until the 
resonance is encountered a second time when the ampli- 
tude decreases. If then > a major part of the 
wave function does not "jump over" the avoided crossing 
back to the initial state, but rather stays in the continu- 
ously connected quasienergy states. This implies that a 
major fraction of the condensate atoms is excited at the 
end of the pulse, escaping out of the lattice. When such 
an experiment is performed repeatedly with fixed rise 
time T p (1) /2 while varying the switch-off duration Tp 2 -*/2, 

one should observe survival probabilities which drop ex- 

(2) 

ponentially with increasing Tp , allowing one to extract 
the quasienergy separation 5s at the avoided crossing 
from the drop rate by means of the known Landau- Zcncr 
formula for quasienergy states [27| ■ 

The results of a series of solutions to the Schrodinger 
equation corresponding to this scenario are plotted in 
Fig. [5l From the slope of the numerical data we de- 
duce a quasienergy gap Se/(hu>) = 0.0099, in agreement 
with the value 0.01 read off from Fig. |4j These find- 
ings clearly underline that this quasienergy gap Se <C huj 
is the actually relevant energy scale for the multiphoton 
transitions under scrutiny here, not the band separation 
Ai = 5.05 hw indicated in Fig. [1] 

There are further options offered by BECs in driven 
optical lattices which have no match in laser-based 
multiphoton studies. For example, one can switch off 
the driving force abruptly at any moment, and analyze 
the state of the wave packet at that particular instant. 
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FIG. 5: Survival probability of atoms in the lowest band after 
asymmetric pulses with fast switch-on time Tp 1 ^ /2 — 5T, and 
with varying switch-off durations Tp 2 '/2. The initial packets 
were centered around (&)(0)/&l = 0.8, so that the dynamics 
are determined by the spectrum shown in Fig. [4] if max = 1.2 
for all pulses. 




T out /T 

FIG. 6: (color online) Pulse tracking realized by switching 
off the squared-sine envelope ((3| with T p = 60 T abruptly at 
t — T ou t, and recording the escape probability at this moment. 
The onset of interband transitions then reveals the presence 
of a multiphoton resonance at the instantaneous amplitude 
reached at t = T out . 



any of the resonances observed in Fig. [3] However, for 
-Kmax = 2.9 interband transitions set in at T out ~ 25 T, 
that is, when the instantaneous amplitude reaches 
K = K max s(T out ) ~ 2.7; this is due to a tiny resonance 
not resolved in Fig. [3] When K niax — 4.0, pronounced 
interband transitions occur at T out ~ 20 T, correspond- 
ing to the instantaneous amplitude K ss 3: This is 
already in the regime of influence of the first of the two 
"large" resonances seen in Fig. [3] Thus, this "pulse 
tracking" strategy allows one to experimentally detect 
multiphoton resonances, that is, to find those values 
of the driving amplitude for which strong resonant 
interband transitions occur. 



IV. SUMMARY AND OUTLOOK 

To summarize, we have argued that BECs in shallow 
optical lattices exposed to ac forcing with frequencies of 
some 100 Hz and smooth envelopes can be employed for 
mimicking multiphoton processes. The enormous degree 
of controllability realizable with such set-ups enables one 
to obtain information not reachable with laser-irradiated 
crystalline solids Q; in particular, we have suggested 
the use of asymmetric pulses for performing avoided- 
quasienergy-crossing spectroscopy. Moreover, we have 
shown how pulse tracking by abruptly switching off the 
driving amplitude allows one to monitor the dynamics 
at each moment during an individual pulse, and thus to 
locate multiphoton resonances between ac-Stark-shifted 
Bloch bands. Obviously, our approach also lends itself 
to systematic explorations of the effects of pulse shaping. 
Even more interestingly, one can activate interparticle in- 
teractions by suitably tuning the s-wave scattering length 
of the Cs atoms [20] • Hence, the detailed experimental 
investigation of the influence of such interactions on mul- 
tiphoton transitions has come into immediate reach. 



Specifically, we again utilize envelopes of the form (|3|), 
and let the amplitude rapidly drop to zero at some 
instant T out , with < T out < T p /2. We set T p = 60 T 
and take initial wave packets with (fe)(0)/&L = 0, as 
for the previous calculations shown in Fig. [2] and plot 
the escape probabilities versus switch-off time T ou t in 
Fig. [HI For -Kmax = 1.6 the wave packet simply follows 
its quasienergy states adiabatically, not encountering 



Acknowledgments 

We are indebted to E. Haller for detailed discussions, 
and for performing preliminary measurements confirm- 
ing the viability of our suggestions. We also acknowledge 
support from the Deutsche Forschungsgcmcinschaft un- 
der grant No. HO 1771/6. 



[1] G. Mainfray and G. Manus, Rep. Prog. Phys. 54, 1333 
(1991). 

[2] C. Cohen- Tannoudji, J. Dupont-Roc, and G. Gryn- 
berg, Atom-Photon Interactions (Wiley-Interscience, 
New York, 1992). 

[3] N. B. Delone and V. P. Krainov, Multiphoton Processes 
in Atoms (Springer- Verlag, Berlin, 2000). 



[4] S. I. Chu and D. A. Telnov, Phys. Rep. 390, 1 (2004). 
[5] S. Arlinghaus and M. Holthaus, Phys. Rev. A 81, 063612 

(2010) . 

[6] S. Arlinghaus and M. Holthaus, Phys. Rev. A 84, 063617 

(2011) . 

[7] N. Gemelke, E. Sarajlic, Y. Bidel, S. Hong, and S. Chu, 
Phys. Rev. Lett. 95, 170404 (2005). 



5 



[8] H. Lignier, C. Sias, D. Ciampini, Y. Singh, A. Zen- 
esini, O. Morsch, and E. Arimondo, Phys. Rev. Lett. 99, 
220403 (2007). 

[9] A. Eckardt, M. Holthaus, H. Lignier, A. Zenesini, D. 
Ciampini, O. Morsch, and E. Arimondo, Phys. Rev. A 
79, 013611 (2009). 
[10] C. Sias, H. Lignier, Y. P. Singh, A. Zenesini, D. Ciampini, 
O. Morsch, and E. Arimondo, Phys. Rev. Lett. 100, 
040404 (2008). 

[11] A. Zenesini, H. Lignier, D. Ciampini, O. Morsch, and E. 
Arimondo, Phys. Rev. Lett. 102, 100403 (2009). 

[12] T. Salger, S. Kling, T. Hecking, C. Geckeler, L. Morales- 
Molina, and M. Weitz, Science 326, 1241 (2009). 

[13] E. Haller, R. Hart, M. J. Mark, J. G. Danzl, L. Reich- 
sollner, and H.-C. Nagerl, Phys. Rev. Lett. 104, 200403 

(2010) . 

[14] J. Struck, C. Olschlager, R. Le Targat, P. Soltan-Panahi, 

A. Eckardt, M. Lewenstein, P. Windpassinger, and K. 

Sengstock, Science 333, 996 (2011). 
[15] R. Ma, M. E. Tai, P. M. Preiss, W. S. Bakr, J. Simon, 

and M. Greiner, Phys. Rev. Lett. 107, 095301 (2011). 
[16] Y.-A. Chen, S. Nascimbene, M. Aidelsburger, M. Atala, 

S. Trotzky, and I. Bloch, Phys. Rev. Lett. 107, 210405 

(2011) . 



[17] J. Struck, C. Olschlager, M. Weinberg, P. Hauke, J. Si- 

monet, A. Eckardt, M. Lewenstein, K. Sengstock, and P. 

Windpassinger, larXiv:1203.0049l 
[18] E. Arimondo, D. Ciampini, A. Eckardt, M. Holthaus, and 

O. Morsch, larXiv:1203. 12591 
[19] O. Morsch and M. K. Oberthaler, Rev. Mod. Phys. 78, 

179 (2006). 

[20] T. Kohler, K. Goral, and P. S. Julienne, Rev. Mod. Phys. 

78, 1311 (2006). 
[21] F. Bloch, Z. Phys. 52, 555 (1929). 
[22] C. Zener, Proc. R. Soc. Lond. A 145, 523 (1934). 
[23] J. Crank and P. Nicolson, Proc. Camb. Phil. Soc. 43, 50 

(1947). 

[24] S. Arlinghaus and M. Holthaus, Phys. Rev. B 84, 054301 
(2011). 

[25] S. Burger, F. S. Cataliotti, C. Fort, F. Minardi, M. In- 

guscio, M. L. Chiofalo, and M. P. Tosi, Phys. Rev. Lett. 

86, 4447 (2001). 
[26] J. Hecker Denschlag, J. E. Simsarian, H. Haffner, C. 

McKenzie, A. Browaeys, D. Cho, K. Helmerson, S. L. 

Rolston, and W. D. Phillips, J. Phys. B: At. Mol. Opt. 

Phys. 35, 3095 (2002). 
[27] K. Drese and M. Holthaus, Eur. Phys. J. D 5, 119 (1999). 



